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Abstract: In this paper, we mainly study the robust stability of linear continuous systems 
with parameter uncertainties, a more general kind of uncertainties for system matrices is 
\ considered, i.e., entries of system matrices are rational functions of uncertain parameters 
which are varying in intervals, we present a method which can check the robust Hurwitz 
stability of such uncertain systems in finite steps. Examples show the efficiency of our 
approach. 
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Given a continuous linear time-invariant system in the state space model, its Hurwitz 
stability is determined by the distribution of eigenvalues of the system matrix. When 
entries of the system matrix are uncertain, e.g., they are varying in intervals, the robust 
stability of such a system have been studied in a large amount of literatures. First, attempts 
were made to find a Kharitonov-like criterion [1] of the stability of an interval matrix 
which only checks some extreme matrices [2], but the criterion was found to be false [3]. 
Later, necessary and sufficient criterions were proposed for interval matrices with special 
properties (e.g., real symmetric interval matrices [7] or Hermitian interval matrices [8]). 
At the same time, various sufficient criterions were found to check the stability of interval 
matrices [4, 5, 6, 7]. 

In this paper, we study a more general kind of uncertainty of system matrices, i.e., 
entries of system matrices are rational functions of uncertain parameters which are bounded 
by intervals, we will present a complete method which can check the robust stability of 
such systems in finite steps. 
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2 Main Results 



Denote by R the field of real numbers, the system matrix A £ M nxn in x(t) = Ax(t) is 
called Hurwitz stable if all its eigenvalues lie in the open left half complex plane. When 
A is continuously varying in M nxn , i.e., A is in a connected set A C M nxn , we say A is 
robustly Hurwitz stable if each A £ A is Hurwitz stable. 

[14] showed that the system matrix with polytopic uncertainty is robustly Hurwitz sta- 
ble if and only if a Hurwitz stable matrix exists and two forms (i.e., homogenous polyno- 
mials) are positive on the standard simplex. In fact, we could come to a similar conclusion 
for A. Denote the characteristic polynomial of A by 

f A (s) = det(sJ B - A) = s n + o^is" -1 + . . . + a lS + a , (1) 

and the Hurwitz matrix of by A^, which is an n x n matrix defined as 

0\ 



' 



V / 

The successive principal minors of are denoted by A&, k — 1, 2, . . . , n. Then we have 

Theorem 1. Suppose some A £ A is Hurwitz stable, then A is robustly Hurwitz stable if 
and only if 

a > and A n _i > for all A £ A. (2) 

Proof. The proof is exactly the same as that of Theorem 1 in [14]. □ 

In this paper, we are interested in a type of matrix uncertainty in which case the entries 
of the matrix are rational functions of parameters varying in intervals, i.e., 

A(q) = (ajj(q)) nxn £ A, 

where q = (qi, . . . , q m ) T , %(q) are rational functions of q, and qu £ [q k , g^], k = 1, . . . , m. 
We have 

Theorem 2. The robust Hurwitz stability of the matrix set A can be checked in finite 
steps. 

The proof of the above theorem will be given in Section 5. 
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3 Simplicial subdivision of the unit hypercubic 

In our method of checking robust Hurwitz stability of A, we need transform this problem 
to a problem of checking positivity of forms on simplices. Since the uncertain parameters 
are varying in hypercubic, we first introduces the procedure [10] of subdividing the unit 
hypercubic [0, l] m into nonoverlapping simplices in this section. 

Denote by m the set of all ml permutations of {1,2,..., m}. Let 9 = {k\k 2 . . . k m ) 6 
6 m , a set of m+ 1 vertexes {a , . . . , a m } of [0, l] m spanning a simplex can be formed using 
following equations. 

a = 0, (3) 
at = ai_i + e fci , i = 1, 2, . . . , m. (4) 

Denote by Se the simplex spanned by {a , . . . , a m }, i.e., 

m m 

S s = {xGl m :x=^ A^, J]) Af = 1, A* > 0, z = 0, . . . , m}, 

i=0 i=0 

it could be readily shown that such constructed Se has the following equivalent definition 

S e = {(xi,..., x m f e R m : 1 > x kl > x k2 > . . . > x km > 0}. 
According to [10], these simplices have no common interior points with each other, and 

[0, l] m = u s e . 

eee m 

4 Positivity of forms on simplices 

Denote by N the set of all nonnegative integers, let a = (ai,a 2 , . . . ,a m ) G N m , and 
|a| = ai + a 2 + ■ ■ ■ + ot m . For a form of degree d 

f[x\ , X 2l ■ ■ ■ , Xm) ^ ^ C a X^ X 2 ' X m , 
\a\=d 

it is immediate that / is strict positive on the standard (m — l)-simplex S m if all c a are 
positive, where 

m 

S m = ■ ■ -,t m ) : y~]tj = l,tj > 0,i = 1, . . . ,m}. 

8=1 

In fact this condition is not only sufficient, but also necessary in the following sense. 
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Theorem 3 (Polya's Theorem, [15]). If a form f(x±, . . . , x m ) is strict positive on S m , then 
for sufficiently large integer N, all coefficients of 

(xi + . . . + x m ) N f(x u . . . 

are positive. 

[16] gave an explicit bound for N, that is 

did- 1) L 

JV >A_J__ ( (, (5 ) 



where 



■ ot\. • ■ ■ Ot m . I I I i , 

L = max <( cj : \a\ = d 

d\ 



and A is the minimum of / on S m . 

A newly proposed method, i.e., the WDS (i.e., weighted difference substitution) method 
[12], can also be used to check positivity of forms efficiently, we will introduce this method 
below. 

Suppose 6 = (kik 2 . . . k m ) G 6 m , let Pe = (pij) mX m be the permutation matrix corre- 
sponding 9, that is 

1, j = h 

o, j^ki • 



Given T m e R mxm , where 



T 



let 



\ ... 

2 m 



V ... i / 



Ao — PeT m , 



(6) 



and call it the WDS matrix determined by the permutation 9. The variable substitution 
x ^— ylgx corresponding 9 is called a WDS. 

In fact, each variable substitution corresponds an assumption of sizes of x i} x 2 , ■ ■ ■ , x m 
in S m . If for each 9 G m , all coefficients of f(Agx) are positive, then /(x) is positive 
on S m . More generally, if there exists k G N, such that all forms in WDS^(/) have no 
nonnegative coefficients, then /(x) is positive on S m , where 

WDS( fc )(/)= |J •■■ U {f{Ae k ---A e ^)} (7) 
is the fcth WDS set of /(x). In fact, the reverse is also true. 
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Theorem 4 ([12]). If f(x u . . . form of degree d, the magnitudes of its coefficients 

are bounded by M, then f is positive on S m , if and only if there exists k < C p (M,m,d), 
such that each form in WDS^ k \f) has no nonnegative coefficients, where 

C P (M, m, d) = 

Remark: The C P (M, m, d) in (8) provides a theoretical upper bound of the number of 
steps of substitutions required to check positivity of an integral form. In practice, numbers 
of steps used are generally much smaller than this bound. 

If coefficients of the form f(x±, . . . all integers with magnitude bounded by M, 

and / is positive on S m , then an explicit positive lower bound of / on S m exists [17], i.e., 

A> {2MY dn n~ dn+1 - d d- nd \ (9) 

The bound in (9) was shown tight in [17], hence from (5) and (8), we know Polya's Theorem 
has a doubly exponential complexity, while the WDS method only has a power exponential 
complexity, as was shown in [14]. 



In (d^ d (d + i)(«-i)("+2)£) - l n A 
In m — ln(m — 1) 



5 Proofs 

Proof of Theorem 2. Let 

Q = {q : Qk e [q k ,qk\, k = i,---,^}> ( 10 ) 

then from Theorem 1, we know that the robust Hurwitz stability of A is equivalent to the 
positivity of rational functions ao(q) and A n _!(q) on Q, which is further equivalent to the 
positivity of polynomials /i(q) and ^(q) on Q, where /i(q) and ^(q) are multiplications 
of the numerators and the denominators of ao(q) and A ra _i(q) respectively. 

Without loss of generality, we can suppose Q = [0, l] m . Otherwise, we can use transla- 
tions and scale transforms of variables in /i(q), /2(q), and obtain new polynomials which 
are required to be positive on [0, l] m . 

The hypercubic [0, l] m can be divided into ml nonoverlapping simplices according to 
the procedure in Section 3, each simplex corresponds a permutation 9j = (jij 2 • • -j m ), 1 < 
j < ml in m , and can be defined as 



S 0j = {q:l> qjl >...>q jm >O}. (11) 

/i(q) or / 2 (q) may be not homogenous on q, if so, we need to homogenize them, i.e., we 
introduce a new variable qo, and let 

hi{qo,qi,...,q m ) = Qo h{— >•••> — ) 

Qo Qo 

if \ degf/2) r /Ql Qm\ 

h 2 {qo,qi, ...,q m ) = Q Q h{— , ■ • • , — )■ 

Qo Qo 
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It is obvious that /i(q) and /2(q) are positive on Sg j if and only if /ii(q) and h 2 (q) are 
positive on Sg j \ {0}, where 

q = (gb,gi, • • • ,?m) T , 

and 

^ = {(go, qi, ■ ■ ■ , Qm) ■ 1 > 9o > > ■ ■ ■ > q Jm > 0}. (12) 

Denote by the unit vector whose kth component is 1 and other components are all 
0, and S m+ i the (m + l)-dimensional simplex in M. m spanned by {0, ei, . . . , e m }, i.e., 

m+l 

S m +i = {(h, ■ ■ • , Wl) : ^ tj < 1, tj > 0, i= 1, ■ ■ ■ , TO + 1}. 

1=1 

Suppose vertexes except of the simplex Sg j are v j0 , . . . , v jm , and the matrix Vj is defined 

as 

= (Vjo, • • • , Vj m ), 

then through a nonsingular linear substitution of variables in /tj(q), i.e., q ^— V^q, we can 
transform Sg j to S'm+i, and obtain a new form ^y(q) = /ij(V^q). It is immediate that h^q) 
is positive on Sg \ {0} if and only if ^^ (q) is positive on S m+ i \ {0}. Since ^(q) has the 
same positivity on S m+ i \ {0} and 5^+1, we finally come to the following result. 

Lemma 1. The matrix set A is robustly Hurwitz stable if and only if following 2m! con- 
ditions are satisfied: 

> for all qG S m+1 , % = 1,2, j = 1,2, . . . ,m\. (13) 

From Equation (5) and (8) in Section 4, we know that the conditions in Lemma 1 can be 
checked in finite steps. Moreover, if the interval vertexes q , q k , k = 1, . . . , m are all rational 

numbers, then coefficients of hij(q) in Lemma 1 are all integers, and from Equation (9), we 
know that the bounds of steps required to check conditions in Lemma 1 can be explicitly 
expressed in m and the coefficient magnitudes and degrees of /i,y(q). □ 



6 Examples 

Example 1. Consider the uncertain system matrix [18] 



( P2 



P2 



2.025 



P-2 



\ 



Pi 



l+pi 



2.025 



where p\ G [l,2],p 2 £ [0,0.5]. Since each of the four forms obtained in (13) has no 
nonnegative coefficients, this system is robustly stable according to Lemma 1. 



6 



References 



[1] V.L. Kharitonov. Asymptotic stability of an equilibrium position of a family of systems 
of linear differential equations. Differentsial'nye Uravneniya 14 (1978) 2086-2088. 1 

[2] S. Bialas. A Necessary and Sufficient Condition for the Stability of Interval Matrices. 
Int. J. Control, 37 (4) (1983) 717-722. 1 

[3] W. C. Karl, J. P. Greschak, G. C. Verghese. Comments on 'A Necessary and Sufficient 
Condition for the Stability of Interval Matrices'. Int. J. Control, 39 (4) (1984) 849-851. 
1 

[4] M. B. Argoun. On Sufficient Conditions for the Stability of Interval Matrices. Int. J. 
Control, 44 (5) (1986) 1245-1250. 1 

[5] Y. T. Juang, S. L. Tung, T. C. Ho. Sufficient condition for asymptotic stability of 
discrete interval systems. Int. J. Control, 49 (5) (1986) 1799-1803. 1 

[6] M. Fu, B. R. Barmish. Maximal unidirectional perturbation bounds for stability of 
polynomials and matrices. Systems & Control letters, 11 (3) (1988) 173-179. 1 

[7] J. Rohn. Positive Definiteness and Stability of Interval Matrices. SIAM J. Matrix 
Anal. Appl. 1994, Vol. 15, No. 1: 175-184. 1 

[8] D. Hertz. The Extreme Eigenvalues and Stability of Hermitian Interval Matrices. IEEE 
Transactions on Circuits and Systems-I: Fundamental Theory and Applications. 1992, 
Vol. 39, No. 6: 463-466. 1 

[9] F. R. Gantmacher. The Theory of Matrices. Chelsea Publishing Company, 1960. 

[10] H. Zhang, S. Wang. Linearly constrained global optimization via piecewise-linear ap- 
proximation. Journal of Computational and Applied Mathematics, 214 (2008) 111-120. 
3 

[11] X. Hou, S. Xu, J. Shao. Some Geometric Properties of Successive Difference Substi- 
tutions. 

[12] X. Hou, J. Shao. Completeness of the WDS method in Checking Positivity of Integral 
Forms. arXiv:0912.1649vl. 4, 5 

[13] B. R. Barmish, M. Fu, S. Saleh. Stability of a Polytope of Matrices: Counterexamples. 
IEEE Transactions On Automatic Control, 33 (6) (1988) 569-572. 

[14] J. Shao, X. Hou. A Complete Method for Checking Hurwitz Stability of a Polytope 
of Matrices, arXiv: 1001.0304. 2, 5 

[15] G. H. Hardy, J. E. Littlewood, G. Polya. Inequalities. Cambridge University Press, 
1934. 4 



7 



[16] V. Powersa, B. Reznickb. A New Bound for Polya's Theorem with Applications to 
Polynomials Positive on Polyhedra. Journal of Pure and Applied Algebra, 164 (2001) 
221-229. 4 

[17] G. Jeronimoa, D. Perruccia. On the minimum of a positive polynomial over the stan- 
dard simplex. Journal of Symbolic Computation, 45 (2010) 434-442. 5 

[18] V. Balakrishnan, S. Boyd and S. Balemi. Branch and bound algorithm for computing 
the minimum stability degree of parameter-dependent linear systems, International 
Journal of Robust and Nonlinear Control 1(4) (1991): 295-317. 6 



8 



